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r*) , Abstract 

<u 

(-H Current observations of the cosmic microwave background could confirm an in- 

'— ^ crease in the radiation energy density after primordial nucleosynthesis but before 

_^^1 photon decoupling. We show that, if the gravitino problem is solved by a light axino, 

^ dark (decoupled) radiation emerges naturally in this period leading to a new upper 

^\ bound on the reheating temperature Tr < 10^^ GeV. In turn, successful thermal 

^-H leptogenesis might predict such an increase. The Large Hadron Collider could en- 

^\J dorse this opportunity. At the same time, axion and axino can naturally form the 
observed dark matter. 
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1 Introduction 

It is a new opportunity to determine the amount of radiation in the Universe from obser- 
vations of the cosmic microwave background (CMB) alone with precision comparable to 
that from big bang nucleosynthesis (BBN). Recent measurements by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) [1], the Atacama Cosmology Telescope (ACT) [2] 
and the South Pole Telescope (SPT) [3] indicate — statistically not significant — the radi- 
ation energy density at the time of photon decoupling to be higher than inferred from 
primordial nucleosynthesis in standard cosmology making use of the Standard Model of 
particle physics, cf. [4, 5]. This could be taken as another hint for physics beyond the 
two standard models. The Planck satellite, which is already taking data, could turn the 
hint into a discovery. 

We should search for explanations from particle physics for such an increase in ra- 
diation [6, 7], especially, because other explanations are missing, if the current mean 
values are accurate. Since tailoring an otherwise unmotivated solution to this particu- 
lar problem is trivial, the origin of the additional radiation should be motivated from 
other problems of cosmology and the Standard Model. Ideally a consistent cosmology is 
possible at the same time and without further ingredients. In this Letter we propose a 
natural origin of the late emergence of the additional radiation of exactly this kind from 
supersymmetry (SUSY) that itself is well known to be motivated independently. 

Supersymmetric theories predict the gravitino as the superpartner of the graviton. 
From gravity mediated SUSY breaking, its mass 7713/2 is expected to set the scale of the 
soft SUSY breaking masses of the superpartner in the minimal supersymmetric standard 
model (MSSM) as msusy ~ O [1^13/2) ■ Since the gravitino interacts only gravitational, 
its decay [8] or decays into the gravitino [9] generically spoil the success of BBN. The 
reheating temperature after inflation Tr were typically restricted to low values excluding, 
for example, thermal leptogenesis [10], which is a simple and elegant explanation for the 
origin of matter with the requirement Tr > O (lO^*^ GeV). In any case, a consistent 
cosmology has to explain the observed matter composition of the Universe and at the 
same time should enable a solution to the strong CP problem of the Standard Model 
with the most attractive solution still found in the Peccei-Quinn (PQ) mechanism [11]. 
The spontaneous breakdown of the PQ symmetry gives rise to a goldstone boson — called 
axion a — accompanied by another scalar — called saxion — and the axino a carrying the 
fermionic degrees of freedom of the axion supermultiplet that has interactions suppressed 
by the PQ scale /a > 6 x 10*^ GeV. 

It has been recognised that an axino of mass m^ < O (keV) [12] as the lightest su- 
persymmetric particle (LSP) with a gravitino next-to-LSP of mass w.3/2 ~ O (100 GeV) 
might provide a natural solution to the cosmological gravitino problem [13]. Since any 
gravitino abundance decays invisibly into axino-axion pairs, reheating temperatures as 
high as Tr ~ 10^^ GeV were claimed to be possible. Due to their suppressed couplings 
axion and axino are indeed decoupled from the thermal bath at such late times when 
the gravitino decays. The axino is light to avoid overproduction. At the same time, 
the lightest ordinary supersymmetric particle (LOSP) can decay into axino and photon. 
Altogether, decay problems and overproduction constraints are circumvented. Recently, 



an extensive investigation of the phenomenological viability of the considered mass hi- 
erarchy, msusy > ""^3/2 > ""^S) reports on numerous restrictions of the PQ and MSSM 
parameter space [14]. In particular, much tighter constraints on the reheating temper- 
ature dependent on the PQ scale and/or arising from the assumed bino LOSP were 
found. 

In this Letter we point out that the reheating temperature is also constrained from 
the allowed amount of relativistic particles decoupled from the thermal bath. Such 
"dark radiation" arises from the gravitino decays, because the axions and axinos with 
iTT-a, "2-a ^ "^3/2 ^^^ emitted relativistically and thus contribute to the radiation energy 
density. This has not been considered so far. Since the amount of dark radiation 
is constrained by different observations, we find a new upper bound on the reheating 
temperature, i.e., Tr < 10^^ GeV, which is four orders of magnitude tighter than the 
original one. This bound is completely independent of the PQ parameter space and its 
dependence on the MSSM parameter space is quite limited. However, any conclusion in 
this Letter relies on the prerequisite that the gravitino is lighter than the LOSP and the 
axino is the LSP. 

The coincidence of the required temperature for successful standard thermal leptoge- 
nesis and the new upper bound on the reheating temperature opens up the opportunity 
to explain the increase in the radiation energy density by the decay of the gravitino, that 
produces unavoidably dark radiation after BBN. We show that at the same time axion 
and axino can naturally form the observed dark matter. 

2 Observational constraints 

Bounds on the radiation energy density are usually given in terms of the effective number 
of neutrino species A'eff defined by 



Prad = 1 + A^cff^ U^ P7 > (1) 




where the radiation energy density /Jrad is given as a sum of the energy density in 
photons p-y = (7r^/15)r^, the energy density in the neutrinos of the standard model 
with A'^^f^ = 3.046 [16] and Ti^/T^ = (4/11)^/^ and any departure from the standard 
scenario parametrised as N^g = N^^ + AA'eff- Since the radiation energy density affects 
the expansion rate of the Universe, AA'efi can be constrained. 

Big bang nucleosynthesis provides a constraint on A'gff consistent with the standard 
model [17] 

N^fP^ = 2.4 ± 0.4 (68% CL). (2) 

This constraint seems to be consistent with earlier studies. A recent study [18] suggests 
an increased N^g at BBN 

AtBBN ^ 3 g ^ Q_35 ^gg^^ (^L) . (3) 

The large central value arises partly from an adapted short neutron lifetime. Further- 
more, the uncertainty is found to be essentially larger in an independent, simultaneous 



study [19] . However, these uncertainties open up the possibihty of an increased A'efr at 
ah times, i.e., < AA'^^^^ = AA'^^^. This is typicahy achieved by stable additional 
degrees of freedom that have once been in thermal equilibrium or the emergence of dark 
radiation from particle decays before BBN. Except comments on the second possibility, 
this will not be discussed in this Letter. We assume that N^ = N^^ = 3.046 and elab- 
orate on the opportunity that A'efr increased after the time of primordial nucleosynthesis, 
i.e., AArJBN ^ g < aA^cm^^ 

The primordial helium abundance is formed at a cosmic time around 1 s and BBN 
ends after roughly 20 min ~ 1.2 x 10^ s. Some 100 ky later structure formation sets 
in, when the matter energy density becomes equal to the radiation energy density at 
teq ~ 4 X 10^^ s, and — on cosmological timescales — shortly after electrons and protons 
recombine to form neutral hydrogen. The Universe becomes transparent for photons. 
We observe the surface of last scattering today as the cosmic microwave background 
(CMB). The current constraints from observations of the CMB are much stronger than 
previous ones [4] but still significantly weaker than those from BBN. Measurements from 
the WMAP satellite [1] using the first and third acoustic peaks and the ground-based 
ACT [2] using observations of the third through the seventh peaks are complementary, 
since they span a broad range of scales. This allows an estimate from the CMB alone [2] 

(CMB alone - ACT) A^^^^^ = 5.3 ± 1.3 (68% CL). (4) 

Recently the ground-based SPT published [3] 

(CMB alone - SPT) A^^^^^ = 3.85 ± 0.62 (68% CL). (5) 

Combined with measurements of today's Hubble expansion rate Hq and baryon acoustic 
oscillations (BAO) the current constraints are [2] 

(ACT + WMAP + BAG +i7o) A^eff^^ = 4-56 ± 0.75 (68% CL). (6) 

and [3] 

(SPT + WMAP + BAO +//o) N^i^^ = 3.86 ±0A2 (68% CL), (7) 

respectively. These limits are consistent with N^g^ at the 2-cr level and any deviation may 
well be due to systematic errors. However, there is a hint of tension and the central values 
are off the standard value by AA^eff ~ 0.81-2.3.^ More important, the Planck satellite 
will significantly increase the precision of CMB observations to AA^gff — 0.26 [21] or even 
better [22]. Thus Planck could reveal a difference between A^|^(= A^^^^^) and A^^ff^^ at 
the level of 3- to 5-0", if the central values from the current measurements are accurate; 
see also [5]. Improvements in the determination of A'^^'^ became crucial. 

To have an effect on the effective number of neutrino species measured from the CMB, 
it is necessary that the additional radiation is generated before the observable modes of 



^ We do not take into account Lyman-a forest data, which may probe N^s at even later times and 
seems to favour even larger central values [20] disfavouring the standard value at 2-a, but also seems to 
be afflicted with large systematic uncertainties. 



the CMB have reentered the horizon. This requirement constrains the maximum hfetime 
^max ^ 1650 y ~ 5.2 x 10^*^ s [23]. This time is before the time of matter-radiation 
equahty, i.e., t^^^ < teq. 

The late emergence of visible radiation like photons is practically excluded. This 
would not only spoil the success of BBN, but also induce a chemical potential for the 
photons that is bounded by CMB observations. 

3 Emergence of dark radiation 

The observations discussed in the previous section constrain the emergence of dark radi- 
ation. Not to affect BBN and to have an effect on measurements of N^s from the CMB, 
the lifetime of a decaying matter particle is constrained to lie in the range 

tfJ^N ~ 1.2 X 10^ s < r < 5.2 X 10^° s ~ t^^^ (8) 

or, equivalently, the cosmic temperature at the particle decay T^j^ ~ 33 keV > T'^^'^ > 
5 eV ~ T'cMB- Since T™j^g < teq in standard cosmology the Universe is radiation domi- 
nated at these times, so that T oc t~^''^. 

We estimate the gravitino decay width with the following considerations. In the 
gravitino center-of-mass frame its two-body decay width is in general given by T^/2 = 
|A^P|pi|/(87rm3,2) . If both decay products are much lighter than the gravitino, their 
equal with opposite sign 3- momenta, |pi|=|p2|, are well approximated by the leading 
term in |pi| = |(m3/2 — O {mi,m2))- The squared amplitude |A^p is on dimensional 
grounds oc ruy^/M^^, where Mpi = 2.4 x 10^^ GeV denotes the reduced Planck mass. 
We have to sum over the outgoing axino spins and average over the incoming gravitino 
spins, which gives a factor 2/4 = 1/2. The rest of the squared amplitude will be a 
factor 1/6 x (1 it c'{ma/m^/2)'^ ± c"{ma/m3/2)^ =t • • •)) where at the leading order the 
gravitino polarisation sum gives a factor 2/3. We assume that a detailed analysis renders 
the typical additional factor 1/4. The coefficients c',c", . . . are usually 1 < c',c" < 12. 
Altogether, we estimate the gravitino decay width T^/2 — 171^,2/ {W2TrM^^) with an 
expected error of the order m-^/m^j2- As r = T^^, we obtain the gravitino lifetime 

V "1-3/2 / 

or, equivalently, T^fg — 24 eV(m3/2/(10^ GeV))'^'^. We see that the gravitino decay is 
expected to happen within the considered range (8). The mass range for the gravitino 
to decay within (8) is 10^ GeV > ^713/2 ^ 35 GeV. Thus most sparticle spectra expected 
at the Large Hadron Collider (LHC) are allowed from this perspective and as soon as 
SUSY is discovered the gravitino next-to-LSP mass is constrained to a small window. 

Since the axino is the LSP, the gravitino next-to-LSP decays into axino and axion. 
Other gravitino channels emitting gluons or photons are never restrictive. Due to an 
additional PQ vertex suppressed by the PQ scale, the branching ratio of such processes 



is extremely small. We find the emitted energy at least twelve orders of magnitudes 
below current limits. 

Qualitatively, this holds true even if i?-parity were broken as long as the breaking is 
not too large. i2-parity violating decays of the gravitino in two standard model parti- 
cles were still gravitational decays and additionally suppressed by the i?-parity violating 
coupling. Therefore, the viability of the late emergence of dark radiation in this sce- 
nario does not depend on conserved i?-parity. We expect the upper bound on i?-parity 
violating couplings from leptogenesis to be more constraining. 

Since the axion and also the axino are much lighter than the decaying gravitino, both 
are emitted with relativistic momenta, pa — Pa — i^z/2/'^- Therefore, the energy density 
in dark radiation p^r consisting of the axion and the axino is given by the energy density 
of the gravitino at its decay. It is 

PdAx^j^ = Pa\Tlj^ + Pa\Tlj^ = P?./2\tIJ^ = "^3/2^3/^(^3/2) ' (^0) 

where the entropy density s = {2it'^ /A:b)g^sT^ could depend on the number of relativistic 
degrees of freedom in the thermal bath g^, and their corresponding temperatures. How- 
ever, for the times considered in this section g^g = 3.91 is fix. If we would not make use 
of the sudden-decay approximation, the resulting energy density would differ by a factor 
■^/7r/2 ~ 0.89 only [24]. The thermally produced gravitino yield Y„^^ is estimated as [25] 

where, as in the following, m-g = rrigdj,) denotes the gluino mass at /x ~ 90 GeV. Since at 
times after BBN g^.^ is fix, the energy density of dark radiation scales equal to the energy 
density of radiation in the thermal bath, even though dark radiation is not coupled to 
the bath. Thus it is most easy to calculate AA'eff at the gravitino decay. 

Using (10) and (1) the change of the effective number of neutrino species by the 
gravitino decay at any time after BBN is given by 
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AA^eff - 0.6 ( ^:^ ) (— ^) 1:^7^,^^^). (12) 



l O^GeV y / rUg y / Tr 
7713/2 ) VlTeV/ V 10^° GeV 

This is a number of O (1) and could have been a priori anything. This coincidence for 
parameter values motivated by completely disconnected reasons is the key observation 
of the present work. We remind that the gluino is expected to be among the heaviest 
sparticles and the gravitino next-to-LSP with rususy ~ "^3/2 ~ 10^ GeV is not only 
motivated from SUSY breaking, but is part of the solution of the gravitino problem as 
described in the introduction. This opens up the opportunity for thermal leptogenesis 
that requires a reheating temperature sufficiently larger than 2 x 10^ GeV. The crucial 
experimental input to determine the minimal temperature are future measurements of 
the light neutrino masses. 

For the approximations made and fixed sparticle masses as in (12) we could state 
that Tr > 9.24 x 10^*^ GeV is excluded at the 5-a level referring to the bound (6) on the 
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Figure 1: The solid lines represent upper bounds on the reheating temperature Tr using (12) 
and (6)as function of the gravitino mass TO3/2 for three different values of the gluino mass 
m-g = 550 GeV, lO"* GeV and 10^ GeV. We restrict 7713/2 < 2 m-g. The shaded regions show 
the corresponding parameter space consistent with our scenario. Indicated are bounds from late 
and early enough gravitino decay (8). As discussed successful thermal leptogenesis requires Tr 
sufficiently larger than 2 x 10^ GeV. The dotted lines correspond to values of AA'^ff = 0.52. 
This is the possible Planck 2-(T exclusion limit, if the observed central value coincides with the 
standard model expectation. The star corresponds to parameter values as appearing in (12). 



effective number of neutrino species. For a fixed reheating temperature the gravitino 
cannot be much lighter, since its mass has the largest exponent in (12) and the gluino 
mass is bounded from below by experiments. The other way around, for a fixed mass 
ratio a heavier gravitino does not allow for much higher reheating temperatures. Thus 
we find an upper bound, Tr < 10^^ GeV, four orders of magnitude tighter than in the 
original work [13]. Different upper bounds on the reheating temperature depending on 
the gravitino mass are depicted in Fig. 1 for different values of the gluino mass. The upper 
bound does not depend at all on the parameters describing the axion multiplet. Even 
if the axino were allowed to be heavy, such that it would not be emitted relativistically, 
the emitted radiation energy would reduce by a factor of 2 only. This is due to the 
axion that still carries half of the energy. In this sense, the upper bound does not rely 
on the very light axino. A dependence on the MSSM parameters enters through the 
gravitino yield (11) only. Since a certain, finite mass gap between gluino and gravitino 



next-to-LSP is well- motivated, we find this dependence quite limited. 

That AA'eff in (12) is of order one opens up another possibility. An increase in N^f^ as 
discussed in the previous section might be predicted by the scenario of successful thermal 
leptogenesis where a light axino solves the gravitino problem, because the, in general 
unknown, reheating temperature becomes practically fixed. ^ This is especially appealing 
in a time where current (on-going) experiments might measure two of the unknowns 
in (12). As mentioned above the Planck satellite mission might reveal a AA^gff > and 
the LHC as collider experiment might measure the gaugino masses, in particular the 
gluino mass. Furthermore, as soon as SUSY is discovered this sheds light on the mass 
of a gravitino assumed lighter than the lightest superpartner in the MSSM. The collider 
phenomenology is not distinguishable from the axino LSP case without gravitino, which 
could itself be mistaken as the neutralino LSP case. However, in contrast to give an 
example heavy charged LOSPs may leave the detector. This would immediately imply 
that dark matter is not formed by neutralinos. Interestingly, the reheating temperature 
in an axino LSP scenario might be probed at the LHC [27, 15]. 

It is important to know at what time the emitted particles become nonrelativistic. In 
general, particles become nonrelativistic at the temperature when their momenta become 
equal to their mass, i.e., p(T^^) ~ m. In the rest frame of a particle decaying into two 
much lighter particles the decay products carry a momentum p{T'^^'^) ~ m/2, where m 
denotes the mass of the decaying particle. After the decay these momenta decrease due 
to the expansion of the Universe. In our scenario this yields for the axino 

Assuming g^:s{Ty^) = g*s{T^'^) we obtain 

-"— '-V(j^)(,^)'. ,14, 

Since this temperature is much smaller than the temperature of the Universe today 
To ~ 0.237 meV, axinos of mass rria ~ O (keV) were surely still relativistic today and 
thus contribute as radiation to today's energy budget of the Universe. For larger axino 
masses rria ~ 200 keV these axinos would become nonrelativistic today and for even 
larger masses rria > 1 MeV these axinos became nonrelativistic in between. They were 
relativistic at photon decoupling and nonrelativistic today. However, the contribution 
of these axinos as dark matter is suppressed by the small mass ratio raa/m^i2- Since 
the axion is much lighter than the axino, as vfia < 10 meV, the emitted axions are still 
relativistic today for all possible parameter values. 

If we require the gravitino to never dominate the energy density of the Universe, 
we find a constraint on Tr similar to the one found above. This is expected, because 



^ In higher-dimensional theories the gravitino density may be independent of Tr for temperatures 
required by thermal leptogenesis [26] . Larger Tr were allowed while the predicted increase in AA'^cff were 
unchanged. 



the Universe is dominated by radiation around the gravitino decay and the decay emits 
a substantial amount of radiation energy. However, there is no significant amount of 
entropy produced in the gravitino decay even if it dominates at its decay, because the 
decay products do not thermahse. Thus cosmological abundances are not diluted and a 
period of matter domination after nucleosynthesis and before matter-radiation equality 
is not excluded by these considerations. Instead, a long period of matter domination is 
excluded by measurements of the radiation energy density as (4), (5), (6) and (7). 

The situation seems more involved for the time window between T™(^g ~ 5.2 x 
10^*^ s and tcq ~ 4 x 10^^ s. It seems reasonable that the release of a huge amount of 
dark radiation from a dominating matter particle at these times affects the CMB in 
an observable way even though maybe only at lower multipoles. An analysis of this 
situation is beyond the scope of this Letter and proposed as future research. 

4 LOSP decay and dark matter 

In the considered scenario neither the gravitino nor the LOSP'^ can account for the 
observed dark matter, but the axion and the axino might be natural candidates. As 
their production depends on the Peccei-Quinn scale fa, we need to consider constraints 
on it from LOSP decay.^ The LOSP has to decay early enough not to spoil the success 
of BBN. We remind that the constraints on the LOSP are greatly relaxed compared to 
the gravitino LSP case, because it decays much faster, since fafniosp/{]^pi'fn3/2)'^ ^ 1- 

Since the lightest neutralino is likely one of the lightest superparticles, we consider 
the case of the neutralino being the LOSP and thus its decay into axino and photon. 
In general, the decay width depends on the bino component, the phase space, the ax- 
ion model and its implementation [24]. For simplicity we assume T(B — )• a -|- 7) = 
agujm~/(1287r'^/^) with Oem = 1/128. Using F = t~^ we obtain the upper bound 
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where we demand the bino lifetime r^ to be at most 0.01 s, which is supposed to be 
conservative. If also the decay into a and Z boson is kinematically unsuppressed, the 
bound is additionally relaxed by a factor of two. Note that a substantial wino or Higgsino 
component lowers the relic density of the neutralino, such that a much later decay 
becomes allowed with r ~ (10^-10^) s [28]. Thus a mixed bino- wino or bino-Higgsino 
state may allow for fa even larger than 10^^ GeV already for neutralino masses close to 
m3/2 ~ 100 GeV. 

The DFSZ [35] axino couples to Higgsino-Higgs and sfermion-fermion via dimension- 
4 operators. Thus Higgsino and stau can be the LOSP. These channels easily allow for 
fa even larger than 10^^ GeV as soon as they are kinematically open. Estimates for any 

^ The lightest ordinary superparticle (LOSP) is the hghtest superpartner in the MSSM. In this sense, 
gravitino and axino are extraordinary superparticles. 

* In this Letter fa = fpci/N is called Peccei-Quinn scale with /pq denoting the scale of PQ symmetry 
breaking if A'^ denotes the number of different vacua. For the KSVZ (DFSZ) model N = 1(G). 



LOSP candidate can be taken from [29] and references therein, respectively."' For LOSP 
masses as large as the upper bound on the gravitino mass, m^M ^ 10^ GeV as following 
from (9), Peccei-Quinn scales as large as 10^^ GeV become allowed already for a bino 
LOSP and other LOSP candidates might even allow for fa ~ 10^^ GeV. 

Axino and axion are produced in the early universe. The axino is produced in thermal 
scatterings after inflation, while the axion is dominantly produced nonthermally from 
the vacuum misalignment mechanism [30] or possible topological defects [31]. Due to its 
tiny mass the thermally produced axion density is negligible. 

In the KSVZ [32] model with heavy exotic quarks the relic axino energy density in 
units of today's critical energy density ^a = '^^'aXaSo / Po with so/po — 5.4 x 10^ GeV~^ 
has been calculated in [24, 33] to 

Oksvz ni r "'M ( ^^ \ /IQiHJeVy 

^^ ~ °-HrtevHi^^^^G^J V~^^J ■ ^ ^ 

This is already close to its density if the axino enters thermal equilibrium 0-^ ~ 1 x 
(ttIq/I keV). A recent study [34] reports that this abundance can be suppressed, if the 
KSVZ quarks are lighter than the scale of PQ symmetry breaking. Furthermore, the 
axino yield is reported to become independent of the reheating temperature in the DFSZ 
model, 

^dfs. ^ 10-6 ^ jg_) (^^^) ' {+^f-n (17) 

for Tr > 10^ GeV, where we also indicate that the DFSZ axino is in addition produced 
via "freeze-in" from heavy Higgsino and/or Higgs decays [36] . The yield from freeze- in is 
again proportional to /~^ and independent of the reheating temperature, but dependent 
on the Higgsino and Higgs masses. In the following three paragraphs we sketch different 
example scenarios to account for the observed dark matter in the considered setting.^ 

/ Natural cold axion dark matter: The density of cold axions from vacuum misalign- 
ment. 



>2 / NU \i 



is of the order of the observed dark matter density S^dm — 0.21 for an initial mis- 
alignment angle O of order one and fa not too far below 10^^ GeV [37]. As we have 
seen, constraints from BBN on the decays of different LOSP candidates often allow for 
large enough Peccei-Quinn scales with the important exception of a bino-like LOSP with 
"^B ^ 500 GeV. At the same time topological defects do not occur, because the scale 
of PQ symmetry restoration is larger than the reheating temperature. From (16) we see 
that the axino is required to be light, i.e., < O (keV) to make up only a small fraction 
of iloM- We mention that an admixture of such light axinos could be favoured from 
problems of small scale structure formation [38]. The DFSZ axino might be allowed 



^ In [29] the decay S — >■ LOSP+. . . has been considered. Thus LOSP and axino need to be interchanged 
in the process. A recent comprehensive study of stau LOSP constraints is found in [15] . 
® Of course this list does not contain all possibilities to achieve a consistent cosmology. 
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heavier, 7n2 ~ 200 keV - 1 MeV, but this depends on the relative contributions from 
the different production mechanisms. Altogether, for fa ~ lO^'^ GeV the dark matter is 
naturally formed by cold axions and the axino is required to be light. 

// Warm axino dark matter: For smaller fa ~ lO^'' GeV the generic axion density 
becomes negligible and the KSVZ axino is overproduced.^ Its mass would be required to 
be at most of O (eV) to satisfy hot dark matter constraints. This situation is disfavoured, 
because we would lack a natural dark matter candidate and models with such small 
masses seem hard to achieve. The DFSZ axino, however, could constitute warm dark 
matter with ttt-q > O (keV) in this case. Thereby it could even be dominantly produced 
by freeze-in. Altogether, for fa ~ 10^*^ GeV the dark matter could be formed by warm 
DFSZ axinos. Much smaller Peccei-Quinn scales are disfavoured for both axion models, 
because the reheating temperature would need to be lowered, so that standard thermal 
leptogenesis were excluded. 

/// Beyond LHC: From the above discussion we have seen, that PQ scales as large 
as 10^^ GeV with correspondingly heavier sparticles are not forbidden for many LOSP 
candidates. Surely, f^™^*^ needs to be suppressed by a vanishing misalignment angle in 
this case and such heavy sparticles are disfavoured in order to stabilise the Higgs mass. 
The KSVZ as well as the DFSZ axino could form cold dark matter in this scenario 
with a much smaller mass gap to the gravitino. In this case, the axino is considerably 
produced from LOSP decays (20), such that a LOSP with a small freeze-out abundance 
is required, and also the gravitino produces a substantial amount of cold axinos. Anyway, 
it is possible — but not arbitrarily — to have the desired dark radiation in form of axions, 
a consistent cosmology and no signal at colliders, because all detectable particles are 
beyond the discovery range of the LHC. 

Besides the production from the discussed mechanisms, the axino LSP is also pro- 
duced in decays of the LOSP. The temperature at which axinos from LOSP decay become 
nonrelativistic depends on the LOSP lifetime riosp as 

where riosp is in good approximation independent of the axino mass. Constraining the 
LOSP lifetime from BBN to be shorter than a second, riosp < 1 s, we see that it is likely 
that axinos from the LOSP decay become nonrelativistic after CMB and before today. 
The contribution of these axinos to today's critical energy density Q.-^"^ '^^^^ is simply 
given by 

(-)losp decay '^a ^fo 

a ~ ~ ^'losp 

''''losp 



"^ For the KSVZ model symmetry restoration would not be problematic here. For the DFSZ model 
we need to ensure that the PQ symmetry is not restored for Tr < 10^'^ GeV. 
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where the mass hierarchy between LOSP and axino can usuaUy not be reduced much. 
The LOSP energy density from freeze-out prior its decay ^[° is highly model and 
parameter dependent. For weakly interacting LOSPs simple estimates tend to ^[°gp ~ 
^DM '^ 0.2. However, we see that even for large values ^1° ~ 10^ x JIdm, as they might 
occur for a bino LOSP [28], axinos from LOSP decay with m^ ^ ?Tiiosp give a negligible 
contribution to fi after they became nonrelativistic. 

On the other hand, axinos from LOSP decay contribute to the radiation energy 
density, as those from gravitino decay, before they become nonrelativistic. We take into 
account the scaling of the energy density of dark radiation pdr oc gJs {T)T^, while we 
simplify exploiting ^^(riosp) = 9*s{tiosp) and assuming g*s(Tiosp) > 9*s{T) = 3.91. The 
resulting change of the effective number of neutrino species is 

AA^i- -- . , X 10-« «,tp (^) ' (S) (^) * , (21) 

where r-^ denotes the fraction of emitted energy that is carried by the axino. If the 
standard model particle emitted together with the axino is also very light or massless 
as the photon, r^ = 0.5 is a good approximation. With the same reasoning as after (20) 
we see that AA'efr from LOSP decay is negligible. This is to some extent self-consistent 
as larger riosp becomes allowed only for smaller r2[°g . 

Saxion Up to now we have not discussed cosmological constraints from and on the 
saxion. Usually the saxion leads to cosmological problems from too late decay as the 
other particles with suppressed couplings. However, since in the considered scenario the 
saxion is allowed to produce LOSPs at any time before BBN^, constraints from and on 
the thermally produced saxion are absent or moderate in all considered dark matter 
scenarios, with the exception of the requirement of a small initial oscillation amplitude 
<^sax ^ fa after inflation [29]. 

In the absence of the LOSP decay problem, the lower bound on the saxion mass from 
early enough decay becomes [29] 

2 

msa. > 108 GeV {jj^) ' , (22) 

where we conservatively neglect the saxion decay into an axion pair. The bound (22) is 
automatically fulfilled for a saxion with rrigax ~ "T-susy as expected, when rrigusy denotes 
the soft SUSY breaking masses. Note that in the considered case of fa ~ 10^^ GeV we 
also consider large values of rususy ^ 10^ GeV. Thus the bound is fulfilled in this case 
as well. 

The saxion decay into two axions, (psax — ?• a + a, is constraint from the requirement 
not to produce too much dark radiation before BBN. In contrast to the late decaying 
gravitino as considered in this work, the saxion has to decay before BBN even if the 



This is the same in the case of broken _R-parity which is discussed in [29] . 
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production of axion pairs is the dominant decay channel. This is due to its nonnegligible 
branching ratio into a gluon pair^. However, it leads to an increase in A'^eff before BBN, 
such that < AiV^ = AN^^ . First, we consider the constraint on the saxion-axion- 
axion self-coupling x by requiring that such an increase should be smaller than one. Then 
the constraint on the branching ratio of the saxion into two axions Baa reads [29] 



B.<0.4(l + 50,^.r^^^j(,S^^j'(-fj^^^-rj , (23) 

where we have approximated the decay width of the saxion as Tsax — Tflx + Tsax- Thus 
our conclusion should hold qualitatively for any axion model. Then the branching ratio 
reduces to 

R 11^ - Z C241 

Since in our scenario the reheating temperature is supposed to be rather large, for 
fa ~ 10^" GeV the saxion enters thermal equilibrium after inflation. In this case the 
saxion yield Y^ax cannot be smaller than the equilibrium value l^tx ~ 1.21 x 10^^. The 
appearing saxion mass rrisax = 10^ GeV in (23) is chosen to show the worst situation 
in the considered scenario. Like a smaller fa or a larger m-sax! also a larger x leads 
to an earlier decay, which corresponds to a smaller rigax at its decay. Thus there is a 
self-curing effect for large x [29]. The bound (23) represents an implicit equation for the 
self-coupling x. Evaluating it for x with parameters fixed as appearing in (23) it turns 
out that there is, indeed, no constraint on x, because branching ratios are by definition 
smaller than one. 

For larger fa ~ 10^^ GeV and fixed Tj^ the saxion does no longer enter thermal 
equilibrium. If we assume its yield equals the axino yield as in (16), we can insert 

into (23). Evaluating again for x we find, that any constraint on x disappears in the 
same sense as before, if msax ^ 290 GeV. Considering /„ ~ 10^^ GeV and larger m-sax 
we see from inserting (25) into (23), that both parameters, if raised like that, lead to a 
vanishing bound on x. 

Even though AA'efj = 1 from saxion decay is consistent with all existing measure- 
ments, for deviations of the effective number of neutrino species AA'^^ < 1.78 it would 
be an unfortunate case considering the discovery potential of Planck. A 3-0" detection 
would become impossible. From the discussion above we see that AA'efr from saxion de- 
cay can well be small for natural parameter values. For m-sax > 900 GeV it is expected 
small in any case and especially for x <C 1 and a; > 1. How to achieve a desired value 
of it, is described in [29]. One can consider cases with a significant deviation from the 



^ Considering a saxion with rrisax ^ 40 MeV, produced from coherent osciUations and Tr < 10® GeV 
the saxion could actually lead to a AA^eff of O (1) after BBN [7]. 
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standard model expectation during BBN and an even larger deviation at photon decou- 
pling, i.e., < AA'^^'^ < AA'^^^. We do not comment further on this possibility. Note 
that if the DFSZ saxion can decay into a Higgs pair, this decay might become stronger 
than the gluon decay channel. This reduces the amount of emitted axions. 

Altogether, we see that we have to consider the saxion decay into two axions, espe- 
cially for the considered scenario. However, depending on the actually measured value 
from Planck this leads to acceptable or mild constraints on the saxion mass. In gen- 
eral, the emergence of dark radiation before BBN would be constrained by a better 
determination of the effective number of neutrino species during BBN. 

Note that bounds on the initial oscillation amplitude are not affected by the absence 
of the NLSP problem. Practically, we assume a mechanism that sets the initial amplitude 

'Psax ~ Ja- 

5 Conclusions and outlook 

The key observation of the present work is: The increase in the effective number of neu- 
trino species in (12) by the invisible gravitino decay after BBN but before photon decou- 
pling is of O (1) for expected (natural) masses and a reheating temperature motivated 
by successful standard thermal leptogenesis, which is a completely disconnected notion. 
Thereby the presupposed mass hierarchy, ?n-susy > "13/2 > ""^S) ^^^ motivated to solve 
the cosmological gravitino problem by — in some sense — the cosmological axino problem. 
Indeed, we have found a new upper bound on the reheating temperature Tr < 10^^ GeV 
in this scenario. Since the decay is gravitational, it occurs generically late and the bound 
does not at all depend on the Peccei-Quinn parameter space. The only requirement on 
the MSSM parameter space stems from thermal gravitino production. It is a large 
enough but finite mass gap between gluino and gravitino next-to-LSP. 

Since neither LOSP nor gravitino are dark matter candidates in the considered sce- 
nario, it is not granted that there are natural dark matter candidates at all. We have 
found that the axion may naturally form the observed amount of cold dark matter with- 
out any conflict with BBN from late LOSP decay. To this end, we commented on the 
BBN constraints from different LOSP candidates extending previous analyses. Estimates 
can be made for any MSSM sparticle using [29]. We mentioned that an admixture of 
thermally produced, light axinos may be favoured from small scale problems of struc- 
ture formation. It seems that the DFSZ axino can form warm dark matter for smaller 
fa ~ 10^*^ GeV. The identification of consistent cosmologies and their corresponding 
parameter space with the desired increase in radiation calls for more comprehensive 
studies. 

We have checked that the scenario is save against other gravitino decay channels and 
radiation from LOSP decay. Constraints from and on the saxion are absent or moderate 
(with the exception of the bound on its initial oscillation amplitude ^g^x ^ fa), because 
it is allowed to produce LOSPs at any time before BBN as in [29]. The most severe 
lower bound on the saxion mass can arise from its decay into an axion pair. We find 
these bounds in the worst scenario still acceptable. The collider phenomenology is not 
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distinguishable from the axino LSP case without gravitino. 

Even though we have also shown that the LHC cannot rule out the scenario, we 
want to point out its surprisingly high testability considering that gravitino and axino 
are generally elusive particles and the reheating temperature is not an experimentally 
accessible quantity. Combining the new opportunity to measure the effective number 
of neutrino species from the CMB alone and the potential discovery of supersymmetry 
at the LHC, two of the unknowns in (12) can be determined and within the scenario 
the gravitino mass would be constrained to a small window as, therefore, the reheating 
temperature as well. The viability of thermal leptogenesis in this scenario were tested 
and the cosmological gravitino problem would turn out as a fortune when the Planck 
mission indeed discovers an increased radiation energy density. In turn, that discovery 
would provide indirect evidence for this particular scenario. 
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